High sensitivity phonon spectroscopy of Bose-Einstein condensates 
using matter-wave interference 
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We study low momentum excitations of a Bose-Einstein condensate using a novel matter-wave 
interference technique. In time-of-flight expansion images we observe strong matter- wave fringe 
patterns. The fringe contrast is a sensitive spectroscopic probe of in-trap phonons and is explained 
by use of a Bogoliubov excitation projection method applied to the rescaled order parameter of the 
expanding condensate. Gross-Pitaevskii simulations agree with the experimental data and confirm 
the validity of the theoretical interpretation. We show that the high sensitivity of this detection 
scheme gives access to the quantized quasiparticle regime. 



The peculiar nature of a Bose-Einstein condensate 
(BEC) is clearly observed in its low momentum excita- 
tions, revealing the many-body coherent nature of the 
system. Recent experiments have provided observations 
of the Bogoliubov excitation spectrum, superfluidity, and 
suppression of low momentum excitations due to quan- 
tum correlations in the ground state 0, 0, Q • 

In these experiments the condensates were typically 
excited by Bragg pulses |j] and imaged after a free ex- 
pansion. At high momenta the excitations were clearly 
separated from the expanding condensate, and could be 
counted and quantified. Even in the phonon regime 
(wavelength larger than but comparable to the healing 
length) a well defined excitation cloud could still be dis- 
tinguished from the condensate (see Fig. la) and was 
found to be amenable to direct atom counting methods. 

However, at sufficiently low momentum (wavelength 
much larger than the healing length), we observe a new 
regime, in which the excitations and the condensate no 
longer separate, regardless of the duration of the expan- 
sion. The excitations are manifested in a clear density 
modulation of the cloud in the absorption images (see 
Fig. lc), reminiscent of in-situ observations of low energy, 
zero momentum, excitations of the condensate 0, Q , but 
they differ significantly from these since they involve both 
nontrivial time-of flight (TOF) dynamics and a small, but 
finite, excitation momentum. 

In this Letter we analyze the dynamics of phonons in 
such freely expanding condensates using both the Gross- 
Pitaevskii equation (GPE) and a dynamically rescaled 
Bogoliubov theory [8j. We find that, in our elongated 
condensate, axial low-momentum phonons are adiabati- 
cally converted by the (mainly radial) expansion into free 
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FIG. 1: Density n(p,z) of the expanding condensate, ini- 
tially excited with phonons of momentum q. In (a) and (b): 
qa p = 2.03. The density in (a) is obtained by computer- 
ized tomography of TOF column density absorption images, 
n co \(y,z), while (b) corresponds to a full GP simulation. In 
both cases, note the clearly separated excitation cloud on the 
right of the condensate. The same quantities are plotted in (c) 
and (d) for qa p = 0.31. Note the strong density modulation, 
with no significant outcoupled fraction. 



atoms with the same axial momentum. The overlap of 
these moving free particles with the expanding ground 
state results in the axial periodic density modulations 
observed in the TOF images. 

We use the fringe visibility of these density modu- 
lations after TOF, as an extremely sensitive spectro- 
scopic probe of the excitation strength. Whereas the 
sensitivity of atom or momentum counting methods Q 
scales quadratically with the quasiparticle excitation am- 
plitude, we show that the fringe visibility scales linearly 
with these amplitudes. This enhanced sensitivity was 
noted and used to estimate the degree of adiabaticity in 
the loading of a BEC into an optical lattice 0], but not 
as a spectroscopic probe of the excitations themselves. 
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Our nearly pure (> 95%) 87 Rb condensate in the 
\F,mf) = |2, 2) ground state, is formed in an elon- 
gated cylindrically-symmetric harmonic trap with ax- 
ial frequency lu z = 2n x 25 Hz and radial frequency 
ui p = 2tt x 220 Hz. The radial harmonic oscillator length 
is a n = (h/mujp) 1 / 2 = 0.73pm. We use condensates with 
N = 10 5 atoms, having chemical potential p — rjhujp, 
with r\ = 9.1, and an average healing length £ ~ 0.3a p . 
Using Bragg excitation beams detuned 6.5 GHz above 
the 55*1/2, F = 2 — ► 5Ps/2>F' — 3 transition, Bogoli- 
ubov quasiparticles are excited in the condensate along 
the z axis, at a controlled wavenumber q by varying the 
angle between the two beams. At a given q the frequency 
difference between the beams to is controlled via acousto- 
optic modulators. After the Bragg pulse (of duration 
ts = 6 msec) the magnetic trap is rapidly turned off and 
after 38 msec of TOF an on-resonant absorption image 
is taken, with the absorption beam perpendicular to the 
z-axis. 

The dynamics of phonons within an expanding conden- 
sate has recently been studied || and shown that, within 
a short time interval (t ~ ^p 1 ), phonons are converted 
into single particle excitations travelling at an axial ve- 
locity of the order of hq/m. This velocity should be com- 
pared to that of the slowly expanding axial boundary, 
which travels at the asymptotic velocity 7ru> 2 Ztf/2u> p , 
where Ztf = a p {u> p / uo z ){2ri) 1 ^' 2 is the Thomas-Fermi ax- 
ial radius of the condensate. Thus at long expansion 
times the excitations can separate out from the conden- 
sate only if qa p > q c a p = 7r(ui z /ujp)(r]/2) 1 / 2 , while for 
q < q c they remain within the condensate at all times. 
For our experimental parameters one has q c a p = 0.76. 
Experimentally, this estimate is remarkably exact. In- 
deed, we begin to observe a distinct excitation cloud at 
qa p ~ 0.72. However, this transition is not sharp and 
a small fraction still overcomes the condensate bound- 
ary even at lower q. This escaped fraction was used to 
find the resonance in our previous measurement of the 
excitation spectrum 0- 

In the phonon regime, but with q > q c , using comput- 
erized tomography y| (see Fig. la), the freely expanding 
excitation cloud was measured to be carrying almost all 
the momentum and excess energy of the system due to 
excitation. By solving the GPE for the TOF dynamics, 
we now confirm this experimental observation quantita- 
tively for this regime. For instance, the lateral cloud in 
Fig. lb carries about 99% of the momentum and more 
than 95% of the energy of the initial phonons. In Fig. lb 
one also notices that the released cloud is radially wider 
than the condensate, in agreement with the experimen- 
tal observation [j| and has an evident "shell structure". 
Both features are associated with the spatial structure 
of the Bogoliubov amplitudes of the in-trap initial ex- 
citations |3| and the presence of excitations with radial 
nodal lines. At this momentum (qa p = 2.03) counting 
the population of the excitation cloud also gives reason- 
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FIG. 2: Linear density 2n pn(p, z)dp measured from TOF 
image with very low momentum excitation qa p — 0.31 (solid 
line). Note the slight asymmetry to right, which is in agree- 
ment with the GPE simulation (dashed line). 

able agreement with the population expected by simple 
structure factor considerations Q. 

Turning our attention to the very low momentum ex- 
citations with qa p = 0.31, we study the TOF density 
modulations of Figs, lc and Id in detail. Fig. 2 shows 
the linear density of the expanded condensate, defined 
as 2ir J pn(p, z)dp. We note a slight asymmetry to the 
right, which is related to the initial direction of propaga- 
tion of the phonon |To| . 
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FIG. 3: Fourier transform of the column density n co i(y,z) 
for the same image as Fig. lc. (a) Fourier transform, (b) 
Integration over central strip of (a). Note the clear separation 
between the central mode and the sidelobes, despite the low 
momentum of the excitation. 

To quantify the fringes of Fig. lc, we take the Fourier 
transform of the measured density profile, and observe 
clear sidelobes (Fig. 3). We define the fringe visibility as 
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the area of the sidelobes divided by the area of the cen- 
tral lobe. This is a powerful "spatial lock- in" method to 
reduce noise in the measurement. It also agrees strictly 
with the standard definition of contrast. We note the 
ease of the measurement, due to the clear separation be- 
tween the Fourier peaks, despite the very low momenta 
involved. This should be contrasted with the atom count- 
ing methods, where it becomes increasingly difficult to 
discern the excitation cloud at finite expansion times in 
real space. At such low momentum we do not expect 
any significant radial excitations, since the coupling to all 
higher radial modes is predicted to be negligible [Tl|, [13 • 
This is confirmed by the lack of any features in the trans- 
verse direction in Figs, lc and 3a. 

Fig. 4 shows the fringe visibility as a function of the 
Bragg frequency detuning (circles). We find the reso- 
nances to be at ±137 ± 10 Hz, which compares favorably 
with the expected value 138 ± 5 Hz for the average Bo- 
goliubov frequency in local density approximation (LDA) 
|lf)| . The width of the resonances is due to the finite du- 
ration of the Bragg pulse. 

A direct GPE simulation of the experiment (Fig. 4, 
solid line) gives good agreement with the observed data 
[l7j . However, little intuition is gained by this. There- 
fore, we also compare this spectrum to the analytic result 
obtained with a rescaled quasiparticle projection method 
@ (dashed line). 

For this, we consider a uniform gas in a box of vol- 
ume L 3 , with a density no = N/L 3 . The ground 
state order parameter is ^0 = \/no, with chemical 
potential /1 = griQ. The excited states are the well 
known Bogoliubov quasiparticles, having frequency Lu k = 
{k 2 /(2m)[h 2 k 2 /(2m) + 2u]} 1/2 and amplitudes u k ,v k = 
±{[h 2 k 2 /(2m) + 2fj]/(2tiw k ) ± 1/2} 1 / 2 . 

We insert the linear expansion 
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into the GPE. The Bragg pulse is included through the 
potential V — Vb cos(qz-ojt), acting for — ts < t < 0. At 
t = —ts the condensate is assumed to be in the ground 
state, so that c k (— is) = for all k's. The equations for 
the coefficients c k (t) are then solved analytically 
After the Bragg pulse one gets 



c±,(0) = T- 



V bV N(u q + V q ) 

2H(u) =F u q ) 



(2) 



and Cfc(0) = for |fc| ^ q. 

Next we assume that the uniform gas undergoes a fic- 
titious radial expansion which mimics the actual expan- 
sion, by taking a uniform time-dependent density that 
decreases in time as the average density of elongated con- 
densate. For this we choose the density no to be equal 



to the radial average of the density of a cylindrical con- 
densate, no = (2/5)n(0), having the same central density 
n(0) of our trapped condensate, and we treat the expan- 
sion by means of the scaling ansatz 0, . 

In this case, the density decreases as n /b 2 (t), where 
b 2 (t) = 1 + uj 2 t 2 is the radial scaling parameter. In the 
presence of quasiparticles with wavevector ±q, one can 
write a rescaled order parameter ^ in the form 



*(z,i) = * +L- 3 / 2 J2 c k (0)u k (t)e ikz +cUO)v* k (t)e- ikz 

k=±q 

(3) 

At t = 0, when G fc (0) = u k , v k (0) = v k and * = ^/n^, 
this expression coincides with QJ. At t > 0, 'J'o re- 
mains stationary while the quantities u and v are time- 
dependent amplitudes that obey rescaled Bogoliubov-like 
equations 
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FIG. 4: Fringe visibility as a function of the Bragg excita- 
tion frequency (points). The error bars represent the uncer- 
tainty due to four measurements. We observe a clear double- 
peaked spectrum, which is finite-time broadened. The peaks 
are found at w = ±137 ± 10, close to the expected Bogoliubov 
frequency (138 ±5 Hz). The solid line is a full GPE simulation 
of the experiment. The dashed line is the result of Eq. (5). 

Further simplification is obtained by assuming adi- 
abatic quasiparticle-to-particle conversion. This corre- 
sponds to assuming that, after an expansion time t ^> 
ujp 1 , the rescaled amplitude u becomes unity, while v 
vanishes ^^|. In terms of ^, this implies 



t) - * + IT 3/2 M0)e^ + c* (OK 



(4) 



Since the expansion is assumed to be radial, the linear 
density obtained by radial integration of the rescaled den- 
sity I 'J 1 2 coincides with that obtained from |^| 2 . The 
same is true for their Fourier components with k — 
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and k — ±q. The sum of the ±q coefficients, di- 
vided by the k = coefficient is exactly the visibility 
7 as defined for the experimental analysis. One finds 
7 = (2/VN)\c q (0) + c* (0)| or, recalling Eq. ©, 



V B 

7 = T 



-1/4 



UJ + UJ q 



(5) 



The dashed line in Fig. 4 corresponds to this expression, 
calculated for the same Vb used in the GPE simulations 
(Vb = O.lfkVp) and with no fitting parameters. The 
agreement with both experiment and GPE simulation in- 
dicates that the essence of the release process of a BEC 
with such very low momentum excitations is indeed cap- 
tured by simply assuming that, due to the reduction of 
the average background density in which the excitations 
live, quasiparticlcs are adiabatically converted into free 
particles that interfere with the ground state. 
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FIG. 5: Fringe visibility after TOF for qa p = 0.31 (cir- 
cles) compared with the standard atom counting technique 
for qa p = 2.03 (boxes), as a function of the number of ex- 
citations. The visibility of the atom counting method grows 
linearly in excitation number with a very weak slope of 1/N, 
whereas the fringe visibility is not negligible (~ 0.4%) even 
for a singly quantized excitation, and grows as the square root 
of the number of excitations. 

We stress that, in this linear response analysis for 
q < q c , the predicted fringe visibility is proportional to 
c q , which is in turn proportional to Vs{u q + v q ). By 
contrast, at higher momentum (q > q c ), the number of 
observed atoms in the excitation cloud is linear in the 



number of excitations, given by 



and quadratic in 



the quasiparticle amplitudes. This is consistent with the 
results of our GPE simulations, where the fringe visibility 
at low q and the counting visibility at high q are found 
to be proportional to Vb and Vj, respectively. 



In our experiments and simulations, we typically ex- 
cite ~ I0 2 quasiparticles. However, by extrapolating the 
expected fringe visibility to lower Bragg intensities, one 
finds that it remains sizeable even in the range of sin- 
gle quantized excitations. We compare (see Fig. 5) the 
fringe visibility at qa p = 0.31 (circles) with the counting 
visibility at q a n = 2.03 (boxes) for our experimental pa- 
rameters |l9j. The fringe visibility at low q is predicted 
to be over two orders of magnitude more sensitive. This 
implies a sensitive heterodyne technique to observe few, 
and possibly single quasi-particle excitations. 

In conclusion, we discuss and implement a novel fringe 
visibility technique to observe the in-situ matter-wave 
interference of low momentum excitations with the con- 
densate. We measure the excitation spectrum with high 
resolution and analyze the possibility of singly quantized 
excitation detection with this technique. 
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between the number of counted atoms in the excitation the total number of atoms N. 

cloud, which is equal to the number of quasiparticles, and 



